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Abstract 

In this paper, colocated MIMO radar waveform design is considered by minimizing the integrated 
sidelobe level (ISL) to obtain beampatterns with lower sidelobe levels than competing methods. First, a 
quadratic programming problem is formulated to design beampatterns by using the minimal-ISL criteria. 
A theorem is derived that provides a closed-form analytical optimal solution that appears to be an 
extension of the Rayleigh quotient minimization for a possibly singular matrix in the quadratic form. 
Such singularities are shown to occur in the problem of interest, but proofs for the optimum solution in 
these singular matrix cases could not be found in the literature. Next, an additional constraint is added 
to obtain beampatterns with desired 3dB beamwidths, resulting in a nonconvex quadratically constrained 
quadratic program which is NP-hard. A semi-definite program and a Gaussian randomized semi-definite 
relaxation are used to determine feasible solutions arbitrary close to the solution to the original problem. 
Theoretical and numerical analyses illustrate the impacts of changing the number of transmitters and 
orthogonal waveforms employed in the designs. Numerical comparisons are conducted to evaluate the 
proposed design approaches. 
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Multiple-input multi-output (MIMO) Radar, transmit beampattem, convex optimization, relaxation. 


I. Introduction 

Multi-input multi-output (MIMO) radars which emit different waveforms from different transmit anten¬ 
nas can provide extra degrees of freedom in waveform design to allow significant coherent gains 0-0 
when the antennas are colocated. Thus, flexible waveform design for transmit heampattem formation is 
of interest for colocated multi-input multi-output (MIMO) radars. Beampattem formation is carried out 
hy designing either a waveform covariance matrix Q, Q or a waveform coefficient matrix |j7|. The goal 
is to focus the transmit power into a range of interesting angles while minimizing the transmit power 
for other angles Q. Various methods can he found in ||T|, Q, ||9[-|[T7| and in the references therein. 
Among the references on heampattem design, the most common strategy is called shape approximation, 
i.e., heampattem design according to a known shape under the criteria of minimum mean square error 
(MMSE) Q, i§, 0, (g, (Tg, i|Tg or minimum difference (MD) Q, To accomplish these 
goals, gradient search algorithms Q. CD. harrier methods CD iterative algorithms ||T4|, ig, and 
convex optimization Q, Q have all been used. Since the sidelobe level is one of the most important 
performance indexes in antenna radiation theory |[T8|-|20|, some minimum sidelobe level design strategies 
have received attention and this is also the topic of this paper. Published methods are the minimum peak 
sidelobe level (PSL) method Q, |12|, |21|, |22|, the discrete prolate spheroidal sequences-based design 
(DPSSD) method |23| and the minimum integrated sidelobe level (ISL) method p^ . The minimum PSL 
method constrains the definition of the sidelobe region as the region outside a band of angles twice as 
wide as the mainlobe, which is a severe restriction that limits flexibility to employ wider mainlobes. 
DPSSD used a maximization of the ratio of energy in mainlobe to the total energy to derive a closed- 
form solution using some approximations. Further, the sidelobe levels obtained by maximizing the ratio 
of energy in mainlobe to total energy can never be smaller than obtained by maximizing the ratio of 
energy in mainlobe to the energy in sidelobes Alternatively, in this paper a closed-form analytical 
solution is given for the criterion of minimum ISL and no approximations are employed. To the best of 
our knowledge, this solution and the proof justifying it, have not appeared in any publications to date. 


Our results appear to be an extension of the Rayleigh quotient minimization |25 Sec. 8.2.3] for a possibly 
singular matrix in the quadratic form. Such singularities are shown to occur in the problem of interest, 
but proofs for the optimum solution in these singular matrix cases could not be found in the literature. 
To augment the closed form results, the minimum ISL criterion is next considered with an additional 
constraint to require a desired 3dB beamwidth, This results in a nonconvex quadratically constrained 
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quadratic program which is NP-hard. A semi-definite program and a Gaussian randomized semi-definite 
relaxation is used to determine feasible solutions arbitrary close to the solution to the original problem. 
Extensive numerical comparisons are conducted to evaluate the proposed design approaches. Additionally, 
a novel theoretical analysis describes the impact of the number of transmitters and the number of transmit 
waveforms on the designed beampattems. 

Notation: Throughout this paper, we use lowercase italic letters to denote scalars, lowercase and 
uppercase letters in bold to denote the vectors and the matrices, respectively. Superscripts of *, T and H 
represent the conjugate, transpose and complex conjugate-transpose (or Hermitian) operators of a matrix, 
respectively, while tr{-} represents the trace of a square matrix and vec{-} represents the vectorization 
operator which creates a column vector from a matrix by stacking its column vectors. || • || denotes the 
Euclidean (or norm of a vector and (g) denotes Kronecker product. The notation I denotes identity 
matrix and 0 denotes all-zero vector or matrix. M, C and Z+ denote the real, complex and positive integer 
spaces, respectively. 


II. Problem formulation 


Consider a colocated uniform linear array (ULA) MIMO radar system equipped with M transmit 
antennas, where the antennas are separated by a half wavelength. Assume each transmitter emits a 
weighted sum of Q independent orthonormal baseband waveforms which are expressed by 0 [t) = 


[(t)i {t) 02 (f) • • • 0Q The weightings are described by a coefficient matrix C = [ci C 2 


'q]’ 


where C G and Cg is used to control the amount of the qth waveform added into a sum at each 

S Q „ „2 

. The total transmit power is fixed at E by setting X/ II 11 = E. Thus the signals 

g=i 

antennas can be described by 


V’(f) = C0) (f). (1) 

After applying a steering vector a{9), the waveforms transmitted to a far-held target at 9 with 9 G 
0 = [—7r/2,7r/2] can be expressed as 

s {t, 9) = {9) C0 (f) (2) 

where a(6') = [1 ••• 

Integrating |s (t, 0) p over a time equal to the support of 0 (t), then the expression for the beampattern 
becomes 0, @ 

P{9) = a.^ {9)CC^a{9). (3) 
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Given the beampattern P{8) expressed in ([^, based on we can write it as 

P{0) = ||C^a(6/)f = ||vec{a^ (0)C}f 
= ||((Iq ® W) vec{C})||^ 

(4) 

= vec^ {C} (Iq (g) (a(0)a^ (0))) vec{C} 

= (Iq (g) A(0)) c = c^A(0)c 

where c = vec{C}, A(0) = a{0)a^ {0) and A{0) = Iq 0 A{0). Note that both A{0) and A{0) are 
Hermitian positive semidefinite matrices for \/0. 

Divide 0, the set of all possible 0, into two disjoint sets called the mainlobe region and 0^/, the 
sidelobe region. To concentrate the transmit power into Q^i as much as possible, we can formulate the 
following minimum integrated sidelobe level (ISL) optimization 


min 

|c||>0 


c^Asic 


s.t. 



(5) 


where 




'0= 


A (0) d9, A^i 4 


'Br, 


A{0)d0. 


( 6 ) 


Problem Q is a quadratic programming problem with a ratio objective (QP-R). When A^i (G 
(pMQxMQ-^ is a positive definite matrix with full rank, then (|^ just can be converted into a Rayleigh 
quotient |[25| Sec. 8.2.3] minimization, and an analytical optimal solution can be obtained |[27|-| 


However, in our case Ami may not be a full rank matrix which means Ami will not be positive definite. 
Here we provide a numerical demonstration of this by considering a case when we vary the width of 
mainlobe from 5° to 160° in Fig. [T] From the figure, we can see that only when the width of Qmi is 
greater than 55° is Ami full rank. Thus, a more general method to obtain the analytical optimal solution 
of Q needs to be considered. However, to the best of our knowledge, we have not found any literature 
giving an analytical solution to Q for Ami not full rank. 

The problem in Q seeks the minimal ISL for the defined &mi and &si- However, in some applications 
a particular 3dB beamwidth is also desired, requiring the addition of some constraints to (|^ to obtain 
the altered optimization by 

c^Asic 


mm 


c^A. 


ml^ 


s.t. i < <1, ve e e.„,. 

2 - c^^A(0o)c “ 


(7) 
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Fig. 1: Demonstration on the rank of A^i vs. angle width of Q^i 


where 0o denotes the maximum power point in Qmi, which is usually chosen to he at the center point 
of the main lohe. The above problem is a nonconvex quadratically constrained quadratic programming 


(QCQP) problem p0| | which is NP-hard. 

In the following section, we will propose the best suitable methods to solve the above two problems, 
respectively. 


III. Waveform desgin for Beampattern formation 
A. Beampattern formation for minimal-ISL only 

Although we have demonstrated that Ami might not be a positive definite matrix, we can confirm fhaf 
Ami + ^si is positive definite as shown in the following Lemma 1. 

Lemma 1: Given = Jq ^ Iq (g) A(0) dO and = Jq ^ Iq ® A{9) dO, where 0^; + Qmi = 
[—7r/2,7r/2], A{6) = a(0)a^(0) and a(0) = [1 Ami + A^; is always 

positive definite. 

Proof of Lemma 1: 
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For Iq (g) A(0) is positive semidefinite, thus given Vc G and ||c|| > 0, we have 

{Asi + Ami) c = f 


Iq ® A{e)dej c 


> c 


H 


= C 


H 


(Iq (g) A(0)) cos(9)d6 c 


Iq (g) / A{0)cos{d)dd c 


( 8 ) 


= (Iq (g) (2Im)) c = 2||c||^ > 0 

Obviously, Agi + Ami is always positive definite, hence, the proof of Lemma 1 is completed. ■ 

Based on Lemma 1, then we can use the following theorem (which can he viewed as an extension 
of the Rayleigh quotient minimization) to get the minimal solution for problem Q(Note that c can 
normalized without loss of any generality for the ratio objective function, thus the constraint ||c|p = E 
in <0 is ignored). 

Theorem I: Define Hermitian matrices A, B G and assume A + B is positive definite while 

A is non-negative definite with rank(A) = ? (0 < <^ < n). Focusing on r G and a minimization 
problem formulated as 


then 


min 

|r||>0 


r^Br 

r-^Ar 


(9) 



\ 

(a.-) 


H 

'B^ (Ai 



( 10 ) 


1) If 0 < ? < n, the minimal solution of Q is given by 

/ 

r = U 

V ■ 

where x^m (£ C‘'^^) denotes the eigenvector corresponding to the minimum eigenvalue of A^ ^ B ^A^ 
B G such that B = Bq - , Ai G is an invertible diagonal matrix ob- 

/ Ai 0 

tained from the eigen-decomposition of A by A = U 

y 0 0 

C^n-‘;)x{n-‘;) ^ submatrices of U^BU given by 

Bq B 2 
B§ Bi 


U^, and Bo G Bi G 


U^BU = 


( 11 ) 


2) If ? = n, the minimal solution of (|^ is given by 

r= fA-V2 


H 


( 12 ) 
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where x^m (£ denotes the eigenvector corresponding to the minimum eigenvalue of A ^/^B(A 

Proof of Theorem 1: 

1) 0 < ? < n 

A is a Hermitian and non-negative matrix with rank(A) = ? hy assumption. Using eigen-decomposition, 
we can factorize r^Ar as 

Ai 0 
0 0 


r^Ar = r^U 


U^r 


= y 


H 


Ai 0 \ (13) 

0 0 

= yi-^Aiyi 

where U is a unitary matrix whose columns are eigenvectors argumented with a set of vectors to form a 
basis and Ai G is an invertible diagonal matrix with positive eigenvalues of A along the diagonal. 
Further, y is defined as 

yi 

y2 


y = U^r = 


Ufr 


Then (|9]) can be formulated as 


r^Br y"U"BUy 

||r||>o r-'^Ar ||yi||>0 yi'^nil-'^yi 


where 11 is invertible and 11 = aJ'^^ for Ai = 


Let X = n^yi. Further simplification of ( [T^ leads to 

yf yf ) u^BU 


r^Br 

mm - = mm 

||r||>o r-'^Ar ||x||>o 


= mm 
||x||>0 


yi 

y2 


x^^x 


x^yf )S 


X 

y2 


x^^x 


(14) 


(15) 


(16) 
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E = 



H 


(17) 






(n“iB2) 


If X is a solution to (|16|, then so is Thus we can focus on the solutions satisfying x^x = 1. 


Then the original prohlem Q is equivalent to 


mm 

||x|P = l 





(18) 


(19) 


The Lagrange function of ( [T^ can he written as 

L(x, y 2 , A) = x^Eqx + yf X + x^S 2 y 2 
+ yfs:iy 2 + A(l- ||x||^). 

To get the stationary points of L(x, y 2 ,A), we use the theorems in pT| , where for a real function 
L of the complex variable vectors ^ and the stationary points are obtained by solving V^L = 0 or 
Vg.L = 0. Here we use the gradients of the conjugates of the vectors to obtain the necessary condition 
for a stationary point as 

Vx*L(x,y 2 , A) = Sqx + S 2 y 2 - Ax = 0 (20) 

Vy*f^(x,y2, A) = 5]^x +Siy2 = 0 (21) 

VAL(x,y 2 , A) = 1 - ||x|p = 0 (22) 

On the other hand, let T = A+B, and left and right multiply the equation with and U, respectively, 
then we have 

u^ru = AU + U^BU 


(23) 
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Since F = A + B is positive definite, then the matrix on the right-hand side (RHS) of (231 is 


positive definite, thus = Bi should also he positive definite. Therefore based on (21 1 we have 
y 2 = —Then combining (201 and (22 1 , all the stationary points should satisfy 


(Hq — ^^ 2 ^) X = Ax 


y2 = ^sfx 


Ixir = 1 


(24) 


To obtain the minimal point of (181, we can find (x,y 2 ,A) among fhe sfafionary points (i.e., the 
constraint of ([24|)) to make the objective function in ([T^ the least by 


mm 

x,y2,A 


s.t. 




(Ho - X = Ax 



(25) 


ya = ^S|^x 


llxir = 1 

Left multiplying the first constraint of (25) by x^ on both sides and inserting the first two constraints 


into the objective function, then (251 can be simplified by 

x^ (Ho-H2Hr's|^)x 


mm 

||x|p=l 


(26) 


Based on the Rayleigh-Ritz theorem p2t Sec. 4.2.2], then the minimal solution of (261 becomes the 


eigenvector corresponding to the minimum eigenvalue of (Hq —), where (Hq —) = 

H 

based on (17 1 . If we assume :x-rnin is the eigenvector, then the mini- 


Bn-BaB^^B? 


a: 


mum value of (0 is achieved by 

r = U 

= U 



V-*--*- J ^min 
H 


-Br'Bf(n-i)"x. 


(27) 


In H 


Xr. 


2)e = n 

When ? = n, A is a positive definite matrix and Ai G Using the same process as for the case 

when 0 < e < n, we obtain that the minimum value of (|9l) is |29| 


= fA-V2 


H 


(28) 
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where x^m is the eigenvector corresponding to the minimal eigenvalue of A ^/^B(A 
Therefore the proof of Theorem 1 is completed. 


B. Beampattern formation for desired beamwidths 

In the following, we will propose different convex methods to solve the problem of (|7]) according to 
the relationship between the number of transmitters M and transmit waveforms Q. 

1) Cases for Q = M: Since Q = M, we can use a Hermitian semi-definite matrix R (G to 

represent CC^ = without loss of any generality. Then (4i can be reformulated as 

M 


p{e) = (Iq ® A{e)) c = cf A(0)c, 

9=1 


M 


(29) 


Thus Q can be reformulated as 


mm 


tr ^ A(0)^CgC^ 
9=1 


tr{A|^R} 


= tr{A(0)R} 


R trjA^'R} 


^ 1 ^ tr{A(0)R} 

2 - tr{A(0o)R} - ’ 

R ^ 0 


(30) 


where A^' = /q A{e)d9, Afj = A{e)dO, Ajni = Iq ® A'f} and A*; = I 


Q 


Afj. Also ^ in 


(301 denotes that R should be a semi-definite matrix. 


IS 


The objective function of (30) is a ratio of two linear functions of R Thus if R is a solution to ( |30l ), so 

_ 


Therefore we can focus on the solution that satisfies tr{AjjJ^R} = 1. Then (30) is equivalent 


to 


min tr 
R 


{a!^r} 

{az^k} 


(31) 


< 1 , ye e Qmi 


s.t. tr <1 Aj(}‘R j = 1 

1 ^ tr{A(6»)R} 

2 - tr{A(0o)R} 

R ^ 0 

Obviously, is a semi-definite programming (SDP) problem for the variable R. Thus we can 
use the convex optimization toolbox CVX |34|, ||35| directly to get the optimal solution. 
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2) Cases for Q M: Since x^Tx = tr{Txx^} for Vx and T, then (|7|) is equivalent to 


Hi 


min tr{As;X} 

c,5C 

s.t. tr{Am;X} = 1 


1 ^ tr{A(.)X} ^ ^ 

2 - tr{A(0o)X} - 


(32) 


Observing 


X = cc^ 

, we find that it is not a convex optimization due to the last nonconvex constraint 


X = cc^ (G Therefore, to solve (32), we consider using semidefinite relaxation (SDR) 

1371. Relaxing the constraint X = cc^ to a convex positive semidefinite constraint X — cc^ ^ 0, we 
can obtain a lower bound on the optimal value of ( [3^ by solving the following convex problem 

min tr {A^/X} 

5C,c 


s.t. tr{Am/X} = l 


1 

- < 
2 “ 


tr 


{a(0)x} 


(33) 


tr 


{a(0o)x} 


< 1, yO G &ml 


X — cc^ ^ 0 

Hence problem ( [3^ is the SDR of (jV]) and can be solved with CVX. Suppose the optimal solution of 
(33) is (X*,c*). This may not be the solution of (jv]). To get an accurate approximate solution for (jv]) 


based on (X*, c*), Gaussian randomization |37| can be used. Generate a sufficient number of samples by 
assuming c is a Gaussian variable with c ~ CA((c*, X* — c*c^), then choose the samples satisfying the 
constraint of (0 to find the best feasible point among them to make the objective function the smallest. 


IV. Numerical Results 

In this section, some theoretical and numerical analyses are performed to show how the variables Q 
(number of waveforms) and M (number of transmitters) impact the designed beampatterns. Next, some 
comparison simulations are conducted to verify the performance our methods. 


A. Beampattern analyses with respect to the number of waveforms Q and transmitters M 

1) Analysis on the impact of different Q on the designed beampatterns: First we give a Lemma which 
will be useful to our analyses. 
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Lemma 2: Given Hermitian non-negative definite matrices A, B G 

r"Br 


such that min = kmin{> 0), then k 

reC"xi ^ 


= min = min where X(g 

reC"xi ^ Ar tr|AA.| 


if there exists kmin{> 0) 
''™) is a 

positive semidefinite Hermitian matrix. 

See the proof in Appendix [a| Given a non-zero vector c_|_ = c^Y' £ where ci G 

(^MQxi we have 



c-p 


Cl -b cf A|^C2 

= cf (^Iq ® A|^^ 

Cl -b ^Iq ® A^^ 

= tr 1 (^Iq (g) A|^ 

) (cicf -bC 2 C^)| 

[c 2 0^(Q_i)]^- Based on (34i, we have 


(34) 




mm 


^Q+i 


^m) c+ 


c;+6C“«5+i)xi (Iq+1 


mm . , 

tr |(Iq 


A^') c+ 

tr { (Iq ® A^) (cicf + C2C §)} 


O (cicf +C2C^)} 


(35) 


> mm —- 

Xo^o tr {(Iq f 


AiyXo} 


Atf) X„} 


where Xq G and the inequality holds because the problem on the RHS is a SDR of the 

problem on the LHS. On the other hand, we have min < min 

^ C+eC“(Q + l)Xl cf!(lQ + l®A™')c+ — C»(lQ(glA™')c 

since setting the extra components in c+ to zero gives the same answer as the problem on the RHS. 


However, no matter what Q and M are, min 
on Lemma 2. Thus 


c"(Iq<S)A^)c 


c+eC“(Q+i)xi )<^+ ceC"'3xi c®(iQ(S)Ajjqc 

values of Q will not have impact on beampatterns obtained by the analytical method proposed in Section 


c c»(lQ(giAS')c 
cf (Iq+i®AJ^)c+ ^ 


= mm 


tr{(Io®A-^)Xo} 
"O tr{(lQ®AS')Xo} 


will always hold based 


c"(lQ(giA^)c 


. Then we can conclude that the 


III-A Fig. 2(a) has verified our conclusion since all the cases with different Q have the same shaped 


beampattern. Now we consider the convex method proposed in Section III-B| To begin our analyses, a 
lemma is provided. 

Lemma 3: Given a real-valued function f{Q) = mintr {(Iq ® B) X} subject to tr {(Iq ® A*) X} < 
bi (g M), i = 1, K, where A*, B (g and X (g are all Hermitian matrice^ then 

f{Q -I- 1) = f{Q) holds for \/Q G Z+. 


'Here we consider bounded set, since the minimization is a SDP with bounded objective, so the minimal value, i.e., the value 
of /, for each Q should always exist. 
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See the proof in Appendix It’s not hard for us to prove that the SDRs of the minimization (|7]) for 
a case with Q + 1 and a case with Q can he expressed in the form of /(Q + 1) and f{Q), respectively, 
as defined in Lemma 3. Therefore, the SDRs for a case with Q + I and a case with Q have the same 
minimal value based on Lemma 3. Since we find the final solutions for different Q based on different 
SDRs but with the same minimal value, we can conclude that the beampattems designed by the convex 
method for different Q will generally be the same except for minor differences that may exist because we 


find the best solution for each case among a Gaussian randomized sample. Fig. 3(a) shows the numerical 
results with different values of Q. We can see that there are only minor differences between the different 
beampattems, which not only shows the correctness of our conclusion but also shows the good precision 
of Gaussian randomization. 

2) Analysis on the impacts of different M on the designed beampattems: Considering a vector c G 
Af^c = [c^ 0]A^^^[c^ 0]^ holds. Thus given c* = [c^ • • • Cq]"^ G we have 

Q 

^ si , 

■'Q 


cf (Iq ® c* = ^ cf A^c, 
9=1 
Q 

= E 

9=1 


-H 


Cf 0]A^+i[cf 0]^ 


(36) 


= c 


H 






where c* G ^(^+^^*5x1 is an extension of c* by padding a zero at the end of each M elements. Therefore, 
the beampattem we obtain at M by Q can be viewed as a special case of M+1 with M zero-components. 
Thus when we increase the number of transmitters, we obtain a beampattem with a lower sidelobe level 
because the case of M + 1 will provide a lower optimal value than the case of M which means a lower 


sidelobe level. Fig. 2(b) shows the simulation results of the beampattems for different values of M. 
From the figure, it can be seen that the sidelobe level of the beampattem decreases when M increases, 
as the analysis indicated. When M is large enough, the mainlobe width of the beampattem is steady at 
the desired width of 30°. However, when M is not, though increasing M can reduce sidelobe levels, 
it will result in a reduction in the 3dB beamwidth. Thus when using the analytical method proposed 


in Section III-A one should choose M carefully in beampattem formation if it is desirable to obtain 


a good trade-off between the 3dB beamwidth and the sidelobe level. Similarly, based on (361 we also 
can convert the case of M for (|7]) into a special case of M -|- 1 though (|7]) has an addtional constraint 
comparing with Q, which again means that increasing M will lower the sidelobe level. Fig. 3(b) has 
verified our conclusion. 
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Fig. 2: Normalized beampatterns designed by the analytical method for a desired 30° width of the main 
lobe, (a) vs. different values of Q-, (b) vs. different values of M. 


B. Comparisons among the different design methods 

1) Comparisons between the two proposed methods: For M = 8 and Q = 3, we consider three types 
of desired beampatterns: the one-mainlobe case, the two-mainlobe case and the three-mainlobe|^ case. In 
each case the desired width of each mainlobe region is 22° while the space between mainlobes is 50° 
for the two multiple mainlobe cases. Fig. shows the numerical results of the beampatterns designed 
using the two proposed methods. From the three subfigures, it is seen that both of the methods can 

^Multi-beam parallel design for complicated multi-mainlobe cases can be employed when using the minimal-ISL only criterion. 


May 6, 2015 


DRAFT 

























































































ACCEPTED BY IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS 


15 


E 

0 ) 


0 ) 

CD 



Angle (degrees) 
(a) 


E 

(0 

0 ) 

CD 



-20 0 20 
Angle (degrees) 
(b) 


Fig. 3: Normalized beampatterns designed by the convex method with a desired 30° beamwidth, (a) vs. 
different values of Q; (b) vs. different values of M. 


provide useful beampatterns. Note that all the beampatterns designed by adding the constraint on the 
3dB beamwidth achieve a 3dB beamwidth close to that desired as expected. The mainlobes obtained 
using the method without the constraint are slightly different as might be expected. 

2) Comparisons with conventional methods: The minimal-ISL only design method is compared with 
the DPSSD method proposed in |231, since both of them can obtain closed-form solutions. Fig. 5(a)| shows 
the numerical results of the two methods (where M = 8, Q = 2, and the width of the defined mainlobe 
region is 22°). We can see that our method obtains a lower sidelobe level, which means the closed-formed 
minimal solution of our method has better performance in terms of lowering sidelobe level. The minimal- 
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ISL design method for the desired 3dB heamwidth is compared with the shape approximation methods 
(SAMs) which also provide the desired 3dB heamwidth. Now consider two specific SAM methods. The 
first one is the SAM proposed in l|^ Sec. III-C] (denoted as ’SAMmmse’)- Though SAMmmse is a 
covariance matrix (R) design method, it can obtain the globally optimal solution by using the criteria 
of minimum MSE. The second method is the SAM proposed in |[7| (denoted as ’SAMMD’)^ which used 
shape approximation by minimizing the maximum difference. Fig. 5(b)| shows the simulation results of a 
desired beampattern with a 40° heamwidth, where M = 8 and Q = 3. From the figure, we can see that 
the beampattern designed by our method has the lowest sidelobe level. 

3) Narrow-beam comparisons with conventional phased-array radar: In some cases it is desired to 
focus energy on a single angle Ot- This can be accomplished by setting A^i = Iq® (a(0t)a^(0t)) in ^ 
and then applying Theorem 1. In Fig. we present some comparisons of results obtained from this slight 
modification of the minimal-ISF only criterion with those obtained from the conventional phased-array 
(CPA) radar p8| , which can obtain narrowly focused beampatterns Q. The three subfigures in Fig. 
show that the beampatterns designed using our method have lower sidelobe levels when compared to 
those obtained by the CPA radar. When we carefully inspect the 3dB beamwidths of the beampatterns, 
we can observe that our method can achieve a mainlobe width which seems to approach to that of the 
conventional one as we increase M. As Fig. illustrates, the differences of the 3dB beamwidths for 
M = 10, M = 20 and M = 100 are less than 0.1°, 0.02° and 0.0012°, respectively. 


V. Conclusion 

Colocated MIMO radar waveform design for transmit beampattern formation by minimizing ISF has 
been considered in this paper. Both analytical and convex design methods using the criteria of minimum 
ISF are proposed, respectively, to obtain beampatterns with lower sidelobe levels with different goals. 
Under the minimum ISF design criterion, both theoretical and numerical analyses have shown that 
the number of waveforms Q doesn’t have impact on the quality of beampatterns while the number 
of transmitters M does. Further, the larger the value of M, the lower the value of the sidelobe level. 
Finally, numerical comparisons have shown our methods can obtain beampatterns with lower sidelobe 
levels than conventional methods. 


Appendix A 
Proof of Femma 2 

Define fc* and X* A 0 such that A:* = eigen-decomposition to factorize X* as X* = 

Without loss of generality, we assume t^iAx^xf} - tr{Ax!x"} for i = 1,..., n and r^m is 
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the minimal solution of min then 


K = 


(37) 


trjBXj ^ Er=itr{Bxixf} 
tr{AX*} ^-^,tr{Ax,xf} 

. tr{Bxjxf} tr{Bxixf} 
i tr{AxjX^} trjAxixf^} 

^ tr{Br™„r^i„} _ 

— . r A W -1 ~ ^min 

Note that X* can he an arhitrary feasible point of min thus min > kmin- On the other 

X;^0 wiAXj- X;^0 trtAXj- 

hand, we have min — ^rnin for ™in is an SDR of min Therefore the proof of 


XbO t'-iAX} 

Lemma 2 is completed. 


rec-xi 


Appendix B 
Proof of Lemma 3 

Write /(Q + 1) as /(Q + 1) = min tr{(lQ+i ® B)X+} s.t. tr{(lQ+i 0 A*) X+} <bi, i = 1,A, 

X-j- ^0 


since if we let X_(_ = 
have 


X 


0 


0 OmxM 


(E C^(‘7+t)^^(Q+t)) then f{Q + 1) = f{Q) holds, hence. 


we 


f{Q + i)<fiQ) 


(38) 


Xi X 


12 


Define X_i_ = | " I > where Xi E and X 2 E , then we 

M2 -2 ' 

can rewrite f{Q + l) as 


where Xo E C^Qx^Q and Xo = 


f{Q + 1) = min tr {(Ig 0 B) Xi} + tr {BX 2 } s.t. 

X_)_ ^0 

tr{(lQ 0 Aj) Xi} + tr{AjX 2 } < k, i = 1,..., A 
= rnin tr | (Ig 0 B) fXi + X 2 ^ | s.t. 

tr| (Ig 0 Aj) (^Xi + X 2 ^ I < 7 = 1,..., A 

X 2 0 


(39) 


0 0 


M(Q-l)xM(Q-l) 


Obviously, the minimization in f[Q + l) (see the RHS of (39l) can be relaxed to the minimization in 


f{Q), thus we have f{Q + 1) > /(Q). Combining (381 then we have f{Q + 1) = f{Q). Therefore, the 
proof of Lemma 3 is completed. ■ 
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Fig. 4: Beampattern comparisons between the two proposed methods, (a) One-mainlobe case, (b) Two- 
mainlobe case, (c) Three-mainlobe case. 
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Fig. 5: Beampattern comparisons between the proposed methods and conventional methods, (a) Compar¬ 
ison between minimal-ISL only design and DPSSD, (b) Comparison between minimum ISL for a 3dB 
beamwidth and SAMs. 
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Fig. 6 : Single-angle focused beampattem comparisons between the minimal-ISL criterion and CPA radar, 
(a) M = 10, (b) M = 20, (c) M = 100. 
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